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ABSTRACT 

We present a microlensing analysis of 61 Chandra observations of 14 quadruply lensed quasars. X- 
ray flux measurements of the individual quasar images give a clean determination of the microlensing 
effects in the lensing galaxy and thus offer a direct assessment of the local fraction of stellar matter 
making up the total integrated mass along the lines of sight through the lensing galaxy. A Bayesian 
analysis of the ensemble of lensing galaxies gives a most likely local stellar fraction of 7%, with the 
other 93% in a smooth, dark matter component, at a mean impact parameter Rc of 6.6 kpc from the 
center of the lensing galaxy. We divide the systems into smaller ensembles based on Rc and find that 
the most likely local stellar fraction varies qualitatively and quantitatively as expected, decreasing as 
a function of Rc. 



1. INTRODUCTION 

Through decades of study on quadruply gravitation- 
ally lensed quasars, it has been well established that sim- 
ple mass models of lensing galaxies — a monopole plus a 
quadrupole — are fairly successful in describing the over- 
all surface density of matter in a lensing galaxy. However, 
these models give no indication of the type of matter 
present, whether it is in a smooth, dark matter com- 
ponent, in a clumpy component like stars or dark mat- 
ter sub-halos, or in s ome combination of the two (e.g., 
Kochanek eraI]|2QQ6 ). 



Our previous study of flux ratio anomalies in 10 sys- 



I'he mass models are smooth by design, but it has 
become clear that some small-scale structure, i.e., some 
dumpiness, must be present in the lensing galaxies. 
Much of the observational evidence for this comes from 
the "flux ratio anomalies" seen in several systems, in 
which the simple mass model correctly predicts the loca- 
tions of the quasar images but fails on the relative fluxes 
of those images (e.g., Metcalf fc Zhao^2QQ2y ^Kochanek , 
^gaal|[2QQ4] (Pooley et al.||2007[ ). The smal 



all-scale 



struc- 

ture further lenses the background quasar with little ef- 
fect on the positions of the images but large effect on 
their brightness. 

Arguments were put forth for two leading candidates 
(stars or dark matter halos) that might constitute this 
small-scale structure. In the case of mi/Mensing, dark 
matter condensations of 10^-10^ Mc?) are responsible 
rtWambsecanss & Paczyhski"1992» 'Witt et al.""1995; ^Ma^l 
& Schneider 1998; Metcalf & Madau 2001; Dalai ^ 
Kochanek 2002; Chiba 2002), whereas in the case of mi- 
crolensing, stars in the lensing galaxy are respo nsible 
dWitt et al.||1995[ [Schechter fc Wambsganssl|2002 ). 
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tems (Paper I; IPooley et al. 12007 ) as well as s tudies of 
individual lenses such as RX J1131 



al.,,2006, .Kochanek et al. 200 7 
Fgm5+080 ( |Fooley eTaT 



1231 (|Blackburne et 
Chartas et al.||2009|) and 



m\ [2m\ [Morgan et al 



2008 ) provide very strong evidence that microlensing is 



the primary cause of the flux ratio anomalies. We have 
shown that the flux ratios are more anomalous in X-rays 
than at optical wavelengths. This is because the optical 
emitting region of the quasar accretion disk is compara- 
ble in angular size to the Einstein radii of the microlens- 
ing stars while the X-ray emitting region is considerably 
smaller. If millilensing were responsible for the anoma- 
lies, there should be no chromatic effect bet ween X-rays 
and optical (contrary to what is seen; see also Blackburne 
[Ct al.,2011J since both regions would be essentially pomt 
sources compared to the Einstein radius of a dark mat- 
ter sub-halo. In addition, we would not expect temporal 
variation of the flux ratios within a human lifetime in the 
case of millilensing, whereas they are naturally expected 
to vary on timescales of months to years in the case of 
microlensing (as is indeed observed). 

Because the X-rays come from a region much smaller 
than the Einstein radii of the microlensing stars, they 
offer a much cleaner signal of microlensing than what is 
available from the optical, which gives a convolution of 
microlensing and the finite size of the optical emitting 
r egion of the quasar 



Schechter & Wambsganss (2002) explored the mi 



crolensing ettects ot ditterent fractional contributions of 
stars and dark matter to the total surface density, espe- 
cially in regard to the probability of strong observable 
microlensing effects on saddle point images. They found 
that the probability of a strong demagnification of a sad- 
dle point image, which is often seen in the observations, 
was relatively low for stellar fractions of 2% and 100% 
but became appreciable for stellar fractions of 5%-25% 
(see, e.g., their Figure 3). They exploited this finding to 
determine the most likely stellar fraction for an ensemble 
of 11 lensing gala xies at the typical impact param eter of 



image formation (Schechter & Wambsganss 



2004). 



They noted, however, that their analysis produced in- 
consistent results unless they assumed that the optical 
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continuum emitting regions had an extended component. 
As we showed in Paper I, this is indeed the case, and it 
comphcates the use of the optical data for dark mat- 
ter determinations. The X-rays, coming from essentiahy 
a point source region as far as the microlenses are con- 
cerned, do not suffer such comphcations and offer a much 
mor e promising avenue. 

In Pooley et all f2009'), we ap plied the technique of 
Schechter ^ Wambsganss (2004) to Chandra X-ray Ob- 
servatory observations of FG 1115+080 and constrained 
the dark matter fraction to ^80%-95% at a character- 
istic distance of ~6 kpc from the center of the lensing 
galaxy. In this work, we extend that analysis to Chandra 
observations of 14 gravitational lenses: HE 0230— 2130 
(1 obs.), MGJ0414+0534 (7 obs.), HE 0435-1223 
(1 obs.), RXJ0911+0551 (2 obs.), SDSS J0924+0219 
(1 obs.), HE 1113-0641 (1 obs.), PC 1115+080 (6 
obs.), RXJ1131-1231 (22 obs.), SDSS 1138+0314 (1 
obs.), H 1413+117 (2 obs.), B 1422+231 (3 obs.), 
WFIJ2026-4536 (1 obs.), WFIJ2033-4723 (1 obs.), 
and Q 2237+0305 (12 obs.). The observations and data 
reduction are described in Sectionj2] Our analysis of the 
X-ray data to obtain fluxes for each of the four images in 
each observation is presented in Sections |3] and |4l The 
Bayesian microlensing analysis is given in Section]5| and 
we discuss the results in Sections HI We summarize our 
findings in Section [7[ 

2. OBSERVATIONS AND DATA REDUCTION 

We utilize publicly available Chandra observations of 
14 X-ray bright quadruply lensed quasars. All data were 
downloaded from the Chandra archive, and reduction 
was performed using the Chandra Interactive Analysis 
of Observations (CIAO) software, version 4.2. All ob- 
servations were taken with the telescope aimpoint on 
the Advanced CCD Imaging Spectrometer (ACIS) S3 
chip. The data were reprocessed using the CALDB 4.3.0 
set of calibration files (gain maps, quantum efficiency, 
quantum efficiency uniformity, effective area) including a 
new bad pixel list made with the acis_run_hotpix tool. 
The reprocessing was done without including the pixel 
randomization that is added during standard process- 
ing. This omission slightly improves the point-spread 
function (PSF). The data were filtered using the stan- 
dard ASCA grades and excluding both bad pixels and 
software- flagged cosmic-ray events. Intervals of strong 
background flaring were searched for, and a few were 
found. In all cases, the flares were mild enough that 
removing the intervals would have decreased the signal 
to noise of the quasar images since it would have re- 
moved substantially more source flux than background 
flux within the small extraction regions. Therefore, we 
did not remove any flaring intervals. The observation 
IDs, dates of observation, and exposure times are given 
in Tabled 

3. ANALYSIS OF X-RAY SPECTRA 

For each observation of each system, we extracted 
events in large regions which enclosed all four images 
of the quasar. We fit the spectra of these events to de- 
termine the total flux Fx, tot detected in each observa- 
tion. Later, as we describe in Section [4j we perform two- 
dimensional image fitting to determine what fractions of 
Fx,tot to assign to individual images. 



The source extraction regions were 4V92 in radius, and 
we extracted background counts from an annulus around 
each system with an inner radius of T!38 and an outer ra- 
dius of 14^76. For each observation, we simultaneously fit^ 
the sou rce and background spectra in Sherpa 4.2 (Free- 



("cstat" in Sherpa) and the IN elder & Meadl ( 1965 ) op 



man et al. 2001) using mod ified ^Cash (1979 ) stati stics 

'Me 

timization method ("simplex" in Sherpa). Both source 
and background were modeled as absorbed, independent 
power laws. The absorption column density was fixed 
at the Galactic value in the direction of the lens based 



on the maps of Dickey fc Lockmanj ( |1990| ). These simple 
spectral models are meant only to reproduce the gross 
X-ray spectral shape for flux estimation and not to test 
for the presence of additional features such as extragalac- 
tic absorption or spectral emission lines; nonetheless, the 
reduced fit statistics indicate more than adequate agree- 
ment between the simple models and the X-ray data. 
Fits were performed over the 0.5-8 keV energy range. 
The best fit power-law index, reduced est at statistic, and 
unabsorbed X-ray flux are reported in Table [l] for each 
observation. 

To calculate the uncertainties in the X-ray fluxes, we 
used the Sherpa tool "sample .energy _flux." This tool 
used 1000 samples of the power-law index and amplitude 
from their normal distributions to calculate 1000 values 
of the 0.5-8 keV flux. The standard deviation of that 
flux sample is reported as the uncertainty on the X-ray 
flux in Tabled 

4. ANALYSIS OF X-RAY IMAGES 



As in our previous work (Pooley et al. 2006, 2007 



2009), we rely on two-dimensional fitting to determme 



the relative X-ray intensities of the four quasar images 
in each system. In this work, we have explored several 
strategies to achieve the best determinations of these rel- 
ative intensities. 

First, we make three sky images of each observation 
at a resolution of 0^.^0492 per pixel using events in the 
0.3-8.0 keV energy range. The first image comprises all 
events in the reprocessed Level 2 event list, and we re- 
fer to this as a "standard" image. The second and third 
images are made after appl ying the subpixe l event reposi- 
tioning (SER) algorithm of Li et al. (|2004l) and either in- 
cluding all split-pixel events or including only the corner- 
split events. We refer to these images as the "SER" im- 
age and the "SER-CO" image, respectively. These SER 
images were made with the aim to improve the determi- 
nation of the X-ray intensities by improving the effective 
Chandra spatial resolution. Their efficacy is discussed 
below. The main idea behind the SER method is that a 
charge cloud split between two or more CCD pixels can 
be better positioned based on the distribution of charge 
among the pixels. The corner-split events provide a more 
precise repositioning of the event but are much fewer in 
number (see Section 4.1). 



The models we use to fit the images consist of a fixed 
background level determined from a large, source- free 
area near the system plus four other components to rep- 
resent the four quasar images. In all cases, the four com- 
ponents are fixed in their relative positions to each other, 
with the absolute position allowed to vary to the best fit 
location. Previously, we used two-dimensional Gaussians 
for these components, with each of the four Gaussians in 
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a fit constrained to have the same full width at half- 
maximum (FWHM) but allowing the specific value of 
the FWHM to vary to the best-fit value. We repeat that 
analysis using Gaussians, and we now also explore other 



choices to model the quasar images. Our aim is to find 
the best representation of the X-ray data that yields the 
smallest uncertainties in the model amplitudes, i.e., the 
individual fluxes of the four quasar images. 



Table 1 

Chandra Observations of 14 Quadruply Lensed Quasars 



System Number System Name 


Date 


Obsid Exp. 


HM Fraction HS Fraction LM Fraction LS Fraction Total Fx/10~^^ 


PL 


Reduc. 






(ks) 


(erg cm~^ ^~^) 


Ind. 


C-stat 



9 

10 



12 
13 
14 



HE 0230-2130 

MG J0414+0534 
MG J0414+0534 
MG J0414+0534 
MG J0414+0534 
MG J0414+0534 
MG J0414+0534 
MG J0414+0534 

HE 0435-1223 



2000 Oct 14 

2000 
2000 
2000 
2000 
2001 
2001 
2002 



Jan 13 
Apr 02 
Aug 16 
Nov 16 
Feb 05 
Nov 09 
Jan 08 



1642 

417 
418 
421 
422 
1628 
3395 
3419 



14.8 

6.6 
7.4 
7.3 
7.5 
9.0 
28.4 
96.7 



2006 Dec 17 7761 10.0 0.167 



RXJ0911+0551 1999 Nov 02 
RXJ0911+0551 2000 Oct 29 

SDSS J0924+0219 2005 Feb 24 

HE 1113-0641 2007 Jan 28 



PG 1115+080 
PG 1115+080 
PG 1115+080 
PG 1115+080 

PG 1115+080 
PG 1115+080 

RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 
RX J1131-1231 

SDSS 1138+0314 

H 1413+117 
H 1413+117 

B 1422+231 
B 1422+231 
B 1422+231 

WFI J2026-4536 



2000 
2000 
2008 
2008 
2009 
2009 

2004 
2006 
2006 
2006 
2006 
2006 
2006 
2007 
2007 
2007 
2007 
2007 
2007 
2007 
2007 
2007 
2008 
2008 
2008 
2008 
2008 
2008 



Jun 02 
Nov 03 
Jan 31 
Nov 02 

Feb 09 
Mar 27 

Apr 12 
Mar 10 
Mar 15 
Apr 15 
Nov 10 
Nov 13 
Dec 17 
Jan 01 
Feb 13 
Feb 18 
Apr 16 
Apr 25 
Jun 04 
Jun 11 
Jul 24 
Jul 30 
Mar 16 
Apr 13 
Apr 23 
Jun 01 
Jul 05 
Nov 11 



419 
1629 

5604 

7760 

363 
1630 
7757 
10730 

10795 
10796 

4814 
6913 
6912 
6914 
6915 
6916 
7786 
7785 
7787 
7788 
7789 
7790 
7791 
7792 
7793 
7794 
9180 
9181 
9237 
9238 
9239 
9240 



2007 Feb 13 7759 18.8 



2000 Apr 19 
2005 Mar 30 

2000 Jun 01 

2001 May 21 
2004 Dec 01 



2007 Ju 



28 



WFIJ2033-4723 2005 Mar 10 



Q 2237+0305 
Q 2237+0305 



2000 Sep 06 

2001 Dec 08 



930 
5645 

367 
1631 
4939 

7758 

5603 

431 
1632 



38.2 
88.9 

28.4 
10.7 
47.7 

10.0 

15.4 

30.3 
9.5 



0.494lg;°/,t 

455 + 0.088 

0.46012:2^8 
0.389tS:°o'56 

411+0.020 



+ 0.029 
0.025 



28.8 
9.8 

17.9 

15.0 

26.5 
9.8 
28.8 
14.6 
14.5 
14.6 

10.0 
4.9 
4.4 
4.9 
4.8 
4.8 
4.9 
4.7 
4.7 
4.4 
4.7 
4.7 
4.7 
4.7 
4.7 
4.7 
14.3 
14.3 
14.3 
14.2 
14.3 
14.3 



0.0861 



0.727lS:lf3 

0.168l°:f3? 



0.585 
0.659 
0.413 
0.348 
0.394 
0.335 

0.634 
0.501 
0.476 
0.418 
0.262 
0.272 
0.279 
0.278 
0.300 
0.272 
0.295 
0.309 
0.347 
0.356 
0.346 
0.339 
0.355 
0.354 
0.353 
0.337 
0.318 
0.327 

0.000 

0.256 



+ 0.043 
0.040 
+ 0.081 
0.072 
+ 0.030 
0.028 
+ 0.076 
0.069 
+ 0.071 
0.064 
+ 0.053 
0.048 

+ 0.035 
0.033 
+ 0.043 
0.039 
+ 0.044 
0.040 
+0.041 
0.037 
+0.014 
0.013 
+ 0.014 
0.014 
+ 0.015 
0.014 
+ 0.017 
0.016 
+ 0.018 
0.017 
+ 0.018 
0.017 
+0.018 
0.017 
+ 0.019 
0.018 
+ 0.020 
0.019 
+ 0.020 
0.019 
+ 0.021 
0.019 
+ 0.017 
0.017 
+ 0.012 
0.012 
+0.010 
0.010 
+ 0.011 
0.011 
+ 0.013 
0.012 
+ 0.011 
0.011 
+0.012 
0.012 



+ 0.048 



+ 0.067 
0.056 



300+0.056 

425+°-°^^ 
414+0-043 

401+° °^^ 

U.^U±_Q Q21 



+ 0.122 
0.121 

+ 0.050 
0.042 

+ 0.039 
0.036 



0.461 

0.192 

0.580 
612+° °^^ 



o.i68i«:sr, 

n 007+0. 079 
U.OOl _o.071 

282+°-°^^ 

PI 97-, +0.067 
-^-0.058 
n 9O/I+0 062 
0.284_Q 

340+0 062 

n QOf^ + 035 

P, O.r>+0.019 



0-216lS°53 0-563 

^+0.122 

-0.082 



0.086 

^00+0.207 



021+° °^^ 

316+0.104 

U.O±U_Q Qgg 



0.112 
0.107 
0.342 
0.401 
0.396 
0.475 



+ 0.023 

0.021 
+ 0.032 

0.028 
+ 0.028 

0.027 
+ 0.078 
-0.070 
+ 0.071 

0.064 
+ 0.066 

0.059 



087+°-°°^ 

U.UO/ _Q QQg 

264+°-°^^ 

U.ZO^_Q Q25 

279+°-°^° 

u.z ' ^-0.027 

350+°-°^'^ 

U.OOU_Q Q33 

624+°-°^^ 

U.OZ^_Q Q27 

627+° °^^ 

U.DZ/ _Q Q27 
U.DUO_Q Q28 

607+° °^^ 
569+° °^° 

U.OUy_Q Q29 

607+°-°^^ 

U.OU/ _Q Q32 

U.OO^ _Q Q3Q 

n 507+0.031 
U.OO/ _Q Q29 
595 + 0.028 
U.OZO_Q Q2g 

C97+0.028 

524+°-°^^ 

U.OZ^_Q Q27 

u.ozy_Q Q24 
476+0.015 

509+°-°^^ 
5-17+0.015 

496+° °^^ 
n 509+0-016 

509+° °^^ 
767+0-173 

'O' _0.163 



+ 0.098 
0.082 



0.477; 

P, 97. +0.052 

p. 01 + 0.022 

U.O±0_Q Q21 

307+0.037 

2Q5 + 0.020 



+0.121 

0.135 

+ 0.069 
0.057 

, + 0.011 
-0.010 

106+°°^^ 

U.±UO_Q Q21 



0.405 
0.322 
0.089 



245+° °^^ 

U.Z40_Q Q35 



0.190 
0.185 
0.196 
0.140 
0.206; 
0.182 
0.197 



0.183l°;o27 0.480; 



+ 0.034 
0.029 
+ 0.034 
-0.029 
+ 0.040 
0.033 
, + 0.025 
0.022 
+ 0.030 
0.026 
, + 0.018 
0.016 
+ 0.009 
0.009 
+ 0.064 
0.056 



171 + 0.037 
-^-0.031 

0.200lS;Sf3 
0.168l°:re 



+ 0.111 
0.122 

+ 0.015 
0.014 
+ 0.022 
0.019 
+ 0.010 
0.009 
+ 0.029 
0.024 
+ 0.020 
0.017 
+ 0.018 
0.015 

, + 0.015 
0.014 
+ 0.018 
0.017 
+ 0.020 
0.018 
+ 0.017 
0.015 
+ 0.007 
-0.006 
+ 0.007 
0.006 
+ 0.007 
-0.007 
+ 0.008 
0.007 
+ 0.008 
-0.007 
+ 0.008 
0.008 
+ 0.008 
0.007 
+ 0.009 
0.009 
+ 0.008 
0.007 
1 + 0.008 
0.007 
+ 0.009 
0.008 
. + 0.008 
0.007 
+ 0.006 
0.006 
+ 0.004 
0.004 
+ 0.005 
0.005 
, + 0.006 
0.005 
+ 0.005 
0.005 
+ 0.005 
-0.005 



o.206iS:°r„ 

P, 90-, +0.045 

U.ZO±_Q Q3g 

247+0-016 
0.267t°:Sfg 
0.287t°:glJ 

104+°-°^^ 
276+°-°^^ 
106+°°^^ 

U.lUO_Q Q-^Q 

100+°°^^ 

U.±UU_Q Qj^g 



0.447 

0.150; 

0.126 
0.120 
0.145 
0.099 
0.098; 

0.220 
0.163 
0.169; 

0.118; 

0.083 
0.083 
0.088; 
O.O8O; 
0.085 
O.O8O; 
0.087; 
0.106 
0.090 
0.083; 
0.103 
0.106 
0.127; 
0.105 
0.106 
0.100 
0.101 
0.101 



0.093 

0.045 
0.111 
0.077 
0.116 
0.065 
0.099 
0.079 



0.170 
0.050 



+0.021 

0.018 
+ 0.014 
0.012 
+ 0.024 
0.020 
+ 0.020 
0.016 
+ 0.022 
0.019 
+ 0.014 
0.012 
+ 0.012 
0.011 
+ 0.005 
0.005 
+ 0.029 
0.025 



v + 0.023 
-0.018 
P) -, 09 + 0.080 
U.l^^_0.054 

0.0841^:^3? 
0.068 



0.154 
0.107 
0.125 
0.106 
0.110 
0.092 

0.059 
0.071 
0.076 
0.113 
0.030 
0.018 
0.028 
0.035 
0.046 
0.041 
0.031 
0.048 
0.039 
0.035 
0.027 
0.026 
0.042 
0.033 
0.024 
0.068 
0.079 
0.063 



+0.101 
0.100 

+ 0.015 
0.014 
+ 0.020 
0.017 
+ 0.010 
0.010 
+ 0.023 
0.020 
+ 0.022 
0.019 
+ 0.017 
0.015 

, + 0.006 
0.006 
+ 0.011 
0.009 
. + 0.011 
0.010 
+ 0.015 
0.014 
, + 0.003 
0.003 
+ 0.003 
0.002 
+ 0.003 
■0.003 
+ 0.004 
0.004 
+ 0.005 
■0.005 
+ 0.005 
0.005 
+ 0.004 
0.004 
+ 0.005 
0.005 
+ 0.005 
0.004 
+ 0.004 
0.004 
.+0.004 
0.004 
+ 0.003 
0.003 
, + 0.003 
0.003 
+ 0.002 
■0.002 
+ 0.002 
0.002 
+ 0.004 
■0.004 
1 + 0.004 
0.004 
+ 0.004 
0.004 



027+°-°^^ 

U.UZ t _o.031 

p, -, ^9 + 0.046 
U.10Z_Q Q3^ 

197+0-040 

u. J-y I _o.034 



0.015 
0.012 
0.017 

0.030 

0.209 

0.225 
0.182 



+ 0.003 
0.003 
+ 0.006 
0.005 
+ 0.003 
0.003 

+ 0.014 
0.012 

+ 0.047 
0.039 

+ 0.018 
0.017 
+ 0.032 
0.028 



2.59 


-|- 





.55 


1.8 


0.82 


7.04 


± 


1 


.22 


1.0 


0.56 


7.38 




1 


.28 


1.0 


0.78 


6.98 


± 


1 


.19 


1.1 


0.55 


8.39 


± 


1 


.38 


1.0 


0.56 


8.08 


-|- 


1 


.21 


1.1 


0.57 


6.74 


-|- 





.63 


1.1 


0.91 


6.79 


-|- 





.33 


1.1 


1.42 


3.76 


± 





.89 


1.9 


0.44 


1.34 


± 





.46 


1.1 


0.74 


1.24 


± 





.99 


1.0 


0.36 


0.47 


± 





.21 


2.2 


0.33 


2.04 


± 





.36 


2.2 


0.47 


4.80 







.56 


1.6 


0.72 


5.16 


± 





.94 


1.7 


0.49 


7.11 


± 





.67 


1.7 


0.71 


3.63 


± 





.86 


1.3 


0.55 


4.05 


± 





.84 


1.3 


0.59 


5.54 




1 


.00 


1.3 


0.67 


20.20 


± 


1 


.36 


1.5 


0.71 


19.56 


± 


1 


.68 


1.7 


0.55 


20.45 


± 


1 


.86 


1.7 


0.56 


17.73 


± 


1 


.77 


1.6 


0.49 


51.37 


± 


2 


.74 


1.8 


0.74 


53.06 


± 


2 


.72 


1.8 


0.78 


49.08 


± 


2 


.69 


1.7 


0.72 


45.30 


± 


2 


.60 


1.8 


0.70 


47.69 


± 


2 


.55 


1.8 


0.65 




_i_ 
It 





.oz 


1 Q 
i.O 


u.oy 


46.85 


± 


2 


.60 


1.8 


0.64 


42.94 




2 


.56 


1.8 


0.68 


44.71 


± 


2 


.45 


1.9 


0.74 


45.17 


± 


2 


.39 


1.9 


0.69 


43.93 


± 


2 


.53 


1.8 


0.63 


57.93 


± 


2 


.92 


1.8 


0.74 


43.47 


-|- 


1 


.48 


1.7 


0.84 


52.51 


± 


1 


.66 


1.7 


0.92 


48.11 


± 


1 


.60 


1.7 


0.84 


36.25 


± 


1 


.38 


1.7 


0.78 


41.61 




1 


.48 


1.7 


0.83 


38.88 


± 


1 


.50 


1.7 


0.78 


1.04 


± 





.50 


0.9 


0.47 


1.22 


± 





.47 


0.4 


0.84 


0.90 


± 





.24 


0.5 


0.84 


10.02 


± 


1 


.42 


1.4 


0.89 


10.56 




2 


.43 


1.5 


0.63 


9.45 


± 


1 


.11 


1.6 


0.91 


5.00 


± 


1 


.05 


2.0 


0.49 


1.33 


± 





.35 


2.1 


0.44 


5.32 


± 





.59 


1.8 


0.72 


4.45 


± 


1 


.24 


1.8 


0.47 
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Table 1 — Continued 



System Number System Name 



Date 



Obsid Exp. HM Fraction HS Fraction LM Fraction LS Fraction 



(ks) 


Q 2237+0305 


2006 


Jan 09 


6831 


7.3 


0.413 


Q 2237+0305 


2006 


May 01 


6832 


7.9 


0.447 


Q 2237+0305 


2006 


May 27 


6833 


8.0 


0.495 


Q 2237+0305 


2006 


Jun 25 


6834 


7.9 


0.511 


Q 2237+0305 


2006 


Jul 21 


6835 


7.9 


0.632 


Q 2237+0305 


2006 


Aug 17 


6836 


7.9 


0.500 


Q 2237+0305 


2006 


Sep 16 


6837 


7.9 


0.505 


Q 2237+0305 


2006 


Oct 09 


6838 


8.0 


0.493 


Q 2237+0305 


2006 


Nov 29 


6839 


7.9 


0.568 


Q 2237+0305 


2007 


Jan 14 


6840 


8.0 


0.547 



Total Fx/10-^^ PL 
(erg cm~^ Ind. 



Reduc. 
C-stat 



+ 0.083 
0.069 
+ 0.059 
0.052 
+ 0.095 
0.079 
+ 0.067 
0.059 
, + 0.082 
0.072 
1 + 0.089 
0.074 
+ 0.086 
0.073 
+ 0.094 
0.078 
+ 0.049 
0.045 
+ 0.055 
0.049 



0.143 
0.136 
0.090 
0.117; 
0.113 
0.137 
0.144 
0.138 
0.116 
0.129 



+ 0.042 
0.034 
+ 0.027 
0.024 
+ 0.032 
■0.026 
+ 0.026 
0.022 
+ 0.025 
0.021 
.+0.037 
0.031 
+ 0.036 
0.030 
+ 0.038 
■0.031 
+ 0.017 
0.015 
+ 0.020 
0.018 



0.331 
0.231 
0.208 
0.216 
0.144; 
0.193 
0.230 
0.204; 
0.198; 
0.171 



+ 0.070 
0.058 
+ 0.036 
0.031 
1 + 0.049 
0.040 
+ 0.034 
0.030 
+ 0.027 
0.023 
+ 0.043 
0.035 
+ 0.046 
0.039 
+ 0.048 
0.039 
+ 0.022 
-0.020 
+ 0.022 
0.020 



0.113 
0.186 
0.206 
0.156 
0.111 
0.171 
0.121 
0.165 
0.118 
0.154 



+ 0.034 
0.027 
+ 0.031 
0.027 
+ 0.050 
0.041 
+ 0.028 
0.024 
+ 0.023 
0.020 
+ 0.039 
0.033 
+ 0.030 
0.025 
+ 0.041 
0.033 
+ 0.016 
0.014 
+ 0.021 
0.019 



3.31 ± 1.15 


1.6 


0.43 


5.51 ± 1.36 


1.7 


0.49 


2.96 ± 1.14 


1.6 


0.37 


6.01 ± 1.54 


1.7 


0.53 


6.62 ± 1.66 


1.5 


0.55 


4.50 ± 1.38 


1.4 


0.47 


3.71 ± 1.55 


1.6 


0.47 


3.37 ± 1.05 


1.7 


0.43 


11.01 ± 1.92 


1.7 


0.54 


7.94 ± 1.65 


1.9 


0.54 



Note. — Columns 5 8 give the fractional contribution of each of 
the HM, HS, LM, and LS images (see Section |5.l| for definitions) to 
the total measured X-ray fiux (given in Column 9). See text for details 
of the spectral and image fitting. 

Our first alternative is ttie P profile, a two-dimensional 
Lorentzian with a varying power law of the form I{r) = 
A{1 + (r/ro)^)~", in which each of the four components 
is constrained to have the same tq and a. Our second 
alternative is a (5 function convolved with an observation- 
specific, ray-traced PSF. Because the Chandra PSF is 
both energy- and position-dependent, a separate PSF is 
constructed for each observation based on the exact off- 
axis location of the four quasar images and the measured 
spectrum of the system in that observation. This infor- 
mation was input to the online Chandra Ray Traceij^to 
produce ray traces of the telescope PSF for each obser- 
vation. These ray traces were projected onto the ACIS 
detector using Marx 4.^ to produce images which were 
then used as the PSF convolution kernel for the 6 func- 
tion fits. 

We therefore have a progression of models in terms 
of shape parameters: a zero-parameter PSF model, a 
one-parameter (FWHM) Gaussian model, and a two- 
parameter (ro and a) /3 model. After the analysis of the 
/3-model fits, we noticed a degeneracy between ro and 
a in all three categories of images (standard, SER, and 




0.4 0.6 0.8 

Best fit Tq in /3 model (") 

Figure 1. Best-fit values of ro and a and la uncertainties. Dark 
blue points are for the /3-model fits to the standard images; medium 
blue is for the fits to the SER image; light blue is for the SER-CO 
image. The black line is a fit to the data taking into account 
uncertainties in both parameters. 
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SER-CO), shown in Figure [T] We fit a straight line to 
the points, taking into account errors in both coordinates 
and found 



a = 3.05(ro/arcsec) + 1.03 



(1) 



We then refit all of the images with a P model where 
a was constrained to follow this relation (hereafter /3c 
because it is constrained), making it essentially a one- 
parameter model. 

We fit each of the four models to each of the three 
classes of images in all 61 observations for a total of 732 
fits. All fits were again performed in Sherpa 4.2 using 
the cstat statistic. We employed a Monte-Carlo based 
optimization method ("moncar" in Sherpa) followed up 
by the simplex method. 

We compare the best-fit position of each of the 12 fits 
for an observation to test for fit fidelity. In general, 
most fits agree in position to better than 0^1, but there 
are some significant outliers, most common in the PSF 
fits, indicating a problem with those fits. Visual inspec- 
tion reveals that the fits with large position discrepan- 
cies found minima in the fit space by zeroing out one or 
more of the four components and shifting the other com- 
ponents to match up with quasar images they are not 
meant to represent. While these were slightly statisti- 
cally better fits, they were not useful representations of 
the data for our purposes. We discard all fits which were 
outliers of 0^25 or more from the rest of the analysis, 
which was 14 out of the 732 fits. 
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Figure 2. Histograms of the fraction of counts from the standard 
imag e used in SER and SER-CO images. The SER algorithm of 
|Li et al.| ( |2004| ) is able to use ~75% of the events on average, while 
the corner-only algorithm uses only ~25% on average. 
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Figure 3. Images of PG 1115+080 ObsID 7757 in the 0.3-8 keV band. The left image is 10^^ on a side and is made at the default resolution 
of 0(^492, which matches the physical pixel size of the ACTS detector. The white box shows the region of the other three panels, which 
are ^" on a side and made at 10 times finer resolution than the default. These images have been smoothed by a Gaussian with a 3-pixel 
FWHM for display purposes only. The color maps are different for the large and small images, but in all cases the intensity scales as the 
square root of the surface brightness from to 2000 counts arcsecond"^. Note that the close pair is severely blended in the standard image 
but has two distinct peaks in the SER and SER-CO images. The loss of counts in the SER-CO algorithm is evident. Black crosses mark 
the best-fit locations of the /5c model to the SER image. See the text for details. 



4.1. Comparison of Standard, SER, and SER-CO 
Results 

As mentioned above, the aim of the SER algorithms 
is to improve the spatial resolution of the X-ray image. 
They do this by using only split-pixel events so there is 
necessarily a loss of signal, and this needs to be weighed 
against the improved resolution. To quantify the signal 
loss, we measure the number of 0.3-8 keV counts in the 
standard image, SER image and SER-CO image in a 6^3 
square region around the lensed quasars. The distribu- 
tions of the fraction of counts in the SER and SER-CO 
images compared to the standard image are shown in 
Figure [2] The SER algorithm is able to utilize about 
75% of the events on average, while the SER-CO image 
can utilize only about 25%. 

The advantage of the SER algorithms is demon- 
strated in Figure 3j which shows a standard image of 
PG 1115+080 made at the default resolution along with 
the standard, SER, and SER-CO images made at 0V0492 
pixel" ^. The blended close pair in the standard image is 
separated into two distinct peaks in the SER image. 

This improvement in resolution can be quantified by 
the width parameters of the Gaussian and P models. Us- 
ing the values from all observations of all systems, we find 
that the best-fit Gaussian FWHM is about 0^06 smaller 
on average in the SER images compared to the standard 
images. It is only about another 0^02 smaller on average 
in the SER-CO images. The best-fit /3-model ro is about 
0'!l2 smaller on average in the SER images compared to 
the standard images, but also only about another 0'!02 
smaller on average in the SER-CO images. The results 
are nearly identical for the /3c model: 0V12 smaller in the 
SER images and only another 0^02 smaller in the SER- 
CO images. This small gain in resolution of ^0^02 of the 
SER-CO images comes at a large price in signal loss. 

We compared the best-fit amplitudes of a specific 
model in the standard, SER, and SER-CO images and 
found reasonably good agreement among all four models. 
There tended to be more outliers in the SER-CO image 
fits, likely due to decreased signal, but most amplitudes 
agreed within la uncertainties. This good agreement is 
reassuring, but the real aim of exploring the SER and 



SER-CO images is the possibility of better constraining 
the amplitudes, i.e., reducing the uncertainty in the best- 
fit model amplitudes. 

In all cases, the model fits to the SER-CO have larger 
amplitude errors on average, again most likely due to the 
large reduction in signal inherent in using that algorithm. 
In the Gaussian, and /3c- model fits, the amplitude er- 
rors in the SER fits and the standard fits are comparable 
(13% for SER versus 12% for standard in the Gaussian 
fits and 15% versus 13% in both the (3- and /3c-model 
fits), whereas the amplitude errors in the SER fits are 
somewhat smaller than those in the standard fits with 
the PSF model (12% versus 17%). 

4.2. Comparison of Gaussian, {3, Ray-traced PSF, and 

Pc Models 

The most important feature of a model is how well it 
represents the data. To explore this for the four mod- 
els, we show histograms of the reduced "cstat" statistic 
in Figure |4j The mean values are almost identical for 
all models, but the Gaussian ones tend to have tails to 
higher values of the reduced statistic than the others. 

Another way to visualize the goodness of the image 
fit is by comparing the radial profiles of the data with 
each of the models. To illustrate this, we choose the ob- 
servation that has the highest number of counts, which 
is ObsID 9181 of RX J1131-1231. The radial profile of 
each quasar image from the center to 0^5 is shown in Fig- 
ure [5j Overlaid are the radial profiles from the best-fit 
models of each type. The Gaussian models are univer- 
sally too squat, and the PSF models are often a poor fit 
in the center. The /3- and /3c-model profiles do the best 
job of matching the data at all radii. Note that these 
are not fits to the radial profiles; rather, they are radial 
profiles of the best-fit models overlaid on radial profiles 
of the data. 

4.3. Discussion of Fits and Choice of Best 
Image/Model Combination 

Considering the fit statistics and radial profiles, the j3- 
and /3c-models appear to be the best choices to represent 
the data. Given the tight correlation seen in Figure [l} it 
is not surprising that both models give nearly identical 
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Figure 4. Histograms of the reduced fit statistic for (top to bot- 
tom) Gaussian model fits, ^^-model fits, PSF fits, and /3c-model fits 
to each of the images (standard, SER, and SER-CO). The arrows 
indicate the mean value of the reduced statistic. For the standard, 
SER, and SER-CO images, respectively, these are 0.21, 0.19, and 
0.10 for the Gaussian fits, 0.19, 0.16, and 0.09 for the /5-model fits, 
0.20, 0.17, and 0.09 for the PSF fits, and 0.19, 0.16, and 0.09 for 
the /3c-model fits. 

results. We have a slight preference for the ^^c- models 
because of the reduction in free parameters. 

Although the SER-CO fits have smaller reduced statis- 
tics on average, they have larger uncertainties and more 
outliers and are therefore not the ideal choice. Between 
the standard and SER images, the uncertainties are com- 
parable, as are the reduced statistics. We favor the SER 
images for the task at hand because we believe that the 
increase in effective spatial resolution will provide higher 
fidelity results for the systems where the separation be- 
tween quasar images is far less than V\ 

The results from the /3c-model fits to the SER images 
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Figure 5. Radial surface brightness profiles (SBPs) of the four 
images (top to bottom: A, B, C, and D) of RX J1131-1231 from 
ObsID 9181. Profiles are made from the (left to right) standard, 
SER, and SER-CO sky images. The data and uncertainties are 
shown as black points with the horizontal bar indicating the bin 
width. Profiles of the four model fits of each image are overlaid. 
In most cases, the (3- and /3c-model profiles are indistinguishable. 
Note that the y-axis is logarithmic. 

are given in Table [l] Each amplitude and uncertainty is 
reported as a fraction of the total, defined as the sum of 
the four amplitudes. In many cases the amplitude errors 
are asymmetric. 

5. DARK MATTER DETERMINATIONS 

Our determination of the fraction of stellar matter that 
makes up the total surface mass density for these systems 
relies on the analysis of microlensin g magnification map s 



and follows the Bayesian methods of [ Pooley et al. 
Our specific method is worked out and discussed 



f 



2QQ9D . 
elow. 



5.1. Microlensing Magnification Maps 

The four images of each quasar are either saddle-points 
or minima of the light travel time surface. We denote the 
higher magnification minimum as the "HM" image and 
the lower magnification minimum as the "LM" image. 
Likewise for the saddle point images, the higher magni- 
fication saddle point is "HS", and the lower magnifica- 
tion saddle point is ^TS. " We have previously modeled 
all of these lens systems ( Pooley et al.||2QQ7 Blackburne 
et al. 2011) to determine the local convergence k and 
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shear 7 for each of these images, which also gives the 
"macrolensing" magnification of each image. These pa- 



rameters fflv en in Table [ 2 | are provi ded by the models 
presented in ( Blackburne"et al.]|2011 ). 



Table 2 

Lensing Galaxy Parameters 



HM HS LM LS 



System 


zi 


Im. 


K 


7 


Magnif. 


Im. 


K 


7 


Magnif. 


Im. 


K 


7 


Magnif. 


Im. 


K 


7 


Magni 


HE 0230-2130 


0.52 


A 


0.472 


0.416 


+9.46 


B 


0.510 


0.587 


-9.57 


C 


0.440 


0.334 


+4.95 


D 


1.070 


0.864 


-1.35 


MG J0414+0534 


0.96 


Ai 


0.481 


0.475 


+22.9 


A2 


0.496 


0.544 


-23.9 


B 


0.478 


0.335 


+6.24 


C 


0.618 


0.684 


-3.11 


HE 0435-1223 


0.46 


C 


0.463 


0.394 


+7.51 


B 


0.520 


0.598 


-7.86 


A 


0.460 


0.390 


+7.17 


D 


0.559 


0.637 


-4.73 


RX J0911+0551 


0.77 


B 


0.575 


0.299 


+ 11.0 


A 


0.633 


0.550 


-5.96 


D 


0.286 


0.055 


+ 1.97 


C 


0.650 


0.568 


-5.00 


SDSS J0924+0219 


0.39 


A 


0.490 


0.440 


+ 15.0 


D 


0.517 


0.557 


-13.0 


B 


0.450 


0.390 


+6.65 


C 


0.546 


0.599 


-6.55 


HE 1113-0641 


0.6"^ 


B 


0.484 


0.450 


+ 15.7 


D 


0.510 


0.548 


-16.6 


A 


0.477 


0.441 


+ 12.6 


C 


0.531 


0.570 


-9.53 


PG 1115+080 


0.31 


Ai 


0.537 


0.405 


+ 19.9 


A2 


0.556 


0.500 


-18.9 


C 


0.472 


0.287 


+5.09 


B 


0.658 


0.643 


-3.37 


RX J1131-1231 


0.30 


B 


0.423 


0.507 


+ 13.2 


A 


0.442 


0.597 


-22.2 


C 


0.422 


0.504 


+ 12.5 


D 


0.834 


0.989 


-1.05 


SDSS 1138+0314 


0.45 


A 


0.465 


0.384 


+7.21 


D 


0.523 


0.614 


-6.69 


C 


0.438 


0.349 


+5.15 


B 


0.578 


0.673 


-3.64 


H 1413+117 


0.8't' 


B 


0.454 


0.359 


+5.91 


A 


0.531 


0.634 


-5.49 


C 


0.441 


0.343 


+5.13 


D 


0.576 


0.680 


-3.54 


B 1422+231 


0.34 


A 


0.371 


0.532 


+8.88 


B 


0.400 


0.666 


-12.0 


C 


0.360 


0.485 


+5.73 


D 


1.530 


1.800 


-0.34 


WFIJ2026-4536 


QA^ 


Ai 


0.499 


0.422 


+ 13.7 


A2 


0.528 


0.557 


-11.4 


B 


0.405 


0.299 


+3.78 


C 


0.579 


0.653 


-4.01 


WFIJ2033-4723 


0.66 


Ai 


0.513 


0.267 


+6.03 


A2 


0.621 


0.638 


-3.80 


B 


0.416 


0.290 


+3.89 


C 


0.650 


0.727 


-2.46 


Q 2237+0305 


0.04 


A 


0.400 


0.400 


+5.00 


D 


0.617 


0.617 


-4.27 


B 


0.385 


0.385 


+4.35 


C 


0.721 


0.721 


-2.26 



References. — jBlackburne et al.| ( |2011| ). 
Note. — ^ Estimated. See the text for details. 

These large-scale lens models can give only the to- 
tal n at the site of each image without regard to the 
form of the matter present. We generate a series of 
12 custom microlensing maps for each image by assum- 
ing that some fraction of z^: is in a clumpy component 
(stars) and the rest is in a smooth component (dark mat- 
ter). We use a logarithmic sequence of stellar fractions 
{Sj): 1.47%, 2.15%, 3.16%, 4.64%, 6.81%, 10%, 14.68%, 
21.5%, 31.62%, 46.4%, 68.13%, and 100%. 

In total, 672 microlensing maps were pro duced us- 
ing the "microlens" ray- tracing code (Wambsganss 1990; 
Wambsganss et al.|1990 Wambsganss 1 999). These mag- 
nitication maps are constructed in the source plane, and 
their centers are referenced to the location of one of the 
quasar images. They show the effects of microlensing 
magnification (due to the sum of all the microimages) for 
a source location anywhere within the map. The mean 
macrolensing magnification, due to the smooth lensing 
potential, has been subtracted off. Each map is 2000 x 
2000 pixels, with an outer scale of 20 rEin and a pixel 
size of 0.01 TEin, where rEin is the Einstein radius of a 
microlensing star of average mass. The stars are drawn 
from a mass function similar to the well-known one of 



Kroupa ( 2001, ). The mass function runs from O.O8M0 to 
l.bMQ with a break at O.5M0 and logarithmic slopes of 
— 1.8 and —2.7 below and above the break, respectively. 
The average mass of a microlensing star is O.247M0, and 
the stellar mass above and below which 50% of the mass 
lies is O.335M0. 

Figure [6] shows portions of each of the four microlensing 
maps (HM, HS, LM, and LS) produced for PG 1115+080 
for stellar fractions of both 10% and 100% to illus- 
trate the differences among the microlensing maps. For 
each map, a histogram of the logarithm of the mag- 
nification values is made and normalized, and this is 
used as the probability distribution for microlensing ef- 
fects P{fiij\Sj) where /j^ij = log;Lo(™cromag^ ,,), i G 
{HM, HS, LM, LS}, and Sj is one of the stellar fractions 
listed above. For convenience, we also define rriij = 
—l-ii^j- These normalized histograms are shown in the 



bottom panels of Figure [6) 

5.2. Bayesian Analysis 

Our goal is to determine the probability of each stellar 
fraction Sj for a lensing galaxy. Our measurements of the 
X-ray fiuxes of the four images divided by their respective 
macrolensing magnifications give four estimates of the 
intrinsic fiux Fx,intr of the quasar. We use conditional 
probability to express P{Sj) as 

P{Sj) = ^P{Sj\A:um.^hs.^lm.^ls)P{^) (2) 

where A* = log;Lo(^x,intr/^norm) ^ud P^i indicates the 
estimate of A' from image i. We choose Fnorm = 
10~^^ erg cm~^ s~^, which has no effect on the analy- 
sis. 

We use Bayes's theorem to express 



P(6'j|AHM, AhS, ^LM, Als) = 



^•-Hb? ^'-ijM? ^"-i^by — 
Z]j ^(^HM, ^HS, ^LM, ^Ls\Sj)Ppr{Sj) 

where Pp^iSj) is the a priori probability of Sj and the 
denominator is a normalization term. We take Ppr{Sj) to 
be uniform and combine it with the denominator as the 
constant A in what follows. We compute it by ensuring 
that EjP{Sj) = l. 

Because the four A'^ are physically distinct, their prob- 
abilities are independent from each other, and we can 
express 



P(^HM, ^HS, ^LM, ^Lsl^j) — ]^^(^i|^i) 



(4) 



Substituting Equations ([3| and Q into Equation we 
arrive at 



P{Sj) = A^l[P{X,\Sj)P{X) 



(5) 



X 



and what remains is to calculate P{P^i\Sj) and P(A'). 

The probability of the intrinsic fiux of a quasar P(A') 
can be determined from the number counts obtained 
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Figure 6. Top: small portions (Vie) of the full microlensing magnification maps for each of the four images of PG 1115+080 for both 
Sj = 10% stars (left) and Sj = 100% stars (right). These 250 x 1000 pixel segments illustrate the microlensing differences due to image 
type and stellar fraction. Middle: normalized histograms of the logarithm of the pixel values (fJ'ij) in each microlensing magnification map. 
Bottom: convolution of those histograms with the probability functions of the X-ray fiux of each image using the data from ObsID 363. 
These give the independent probability distributions for the intrinsic fiux of the quasar, A' = logio(^x,intr/10~"'^^ erg cm~^ s"-*^). Plotted 
in the inset is their product Gj. See the text for details. 



from deep stud ies of the X-ray background. We use 
the results from [Giacconi et al.[ ([2001) that N{> Fx) ^ 
^^-0.85^ but we note this has little impact on the anal- 
ysis. A uniform distribution would produce nearly iden- 
tical results. 

We estimate the intrinsic flux Fx,intr from the mea- 
sured flux (/x,i) of an image and the lensing effects, both 
macrolensing and microlensing. The measured flux is 

fx,^ = i^X,intr X X 10^-' (6) 

where A4i is the macro-magnification of image i. We 



define 

Xi = logio([/x,i/i^norm]/A^i) (7) 

which allows us to write 

A'i = + ruij . (8) 

Because the probability of the sum of two random vari- 
ables is the convolution of their individual probabilities, 
we can express 

P{Xi\Sj) = P{xi+mij\Sj) 

= P{xi) * P{mij\Sj) ^' 
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Figure 7. Probability distributions for the stellar fraction, Sj, at the characteristic radial distance Rc from the center of the lensing 
galaxy for 14 quadruply lensed quasar systems. Those labeled in italics do not have a measured lens redshift zi. 



where P{xi) comes from the uncertainties on the flux 
measurements of the images and the P{mij\Sj) are the 
reverse of P{/jjij\Sj), the normahzed histograms of the 
microlensing maps, as discussed above. 

We do not measure the /x,i directly, though; rather, we 
obtain them by multiplying the total flux of all four im- 
ages (via spectral fitting) and the individual fractions of 
the total (via two-dimensional image fitting). We define 



Again, using the property of the sum of two random vari- 
ables, we express 



Pin + T) 
Pin) * P{T) 



(13) 



and 



so that 



T = logio(-Fx,tot/f„orm) 

Ti = logio(fraCi/7Wi) 
Xi = ri+T 



(10) 

(11) 
(12) 



where the probability distributions are asymmetric 
Gaussians with standard deviations equal to the la un- 
certainties in the image fractions (Table jl]) and T is a 
symmetric Gaussian with a standard deviation equal to 
the uncertainty in /x,tot (Table Introducing notation 



(14) 



Gj and using Equations ^ ana(13), we have 
Gj=YlP{X^\Sj) 

i 

= llP{n)*P{T)*P{miJSj) 
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which can be seen in the bottom panels of Figure |6] using 
values from ObsID 363. 

All of the above has been worked out for a single ob- 
servation of a system, but several systems have been ob- 
served multiple times with Chandra. We combine these 
multiple observations using conditional probability: 



P(5,) = ^P(5,|obsfe)P(obsfe) 



(15) 



where we take P{ohsk) as a weighting factor (normalized 
to unity) that combines two measures of the effectiveness 
of the observation to provide unique and useful informa- 
tion. 

The first ingredient in P(obs/c) concerns the unique- 
ness of the information from the observation. Over time, 
the proper motions of the lensing galaxy and background 
quasar, as well as the internal motions of the microlens- 
ing stars, can be thought of as an effective motio n of the 
source through the field of the microlensing map (| Wyithe| 
^ al. 2000). The more time between observations, the 
gher the chance that the source is in a different enough 
region of the map to be considered an independent sam- 
pling of it. We therefore include a term in P(obs/c) pro- 
portional to how isolated in time the observation is, de- 
fined as the sum of the intervals between the observation 
and all other observations. 

The second ingredient in P{ohsk) is based on the qual- 
ity of the information that the observation provides. Ob- 
servations which yield tight constraints on the individual 
fractions and the total flux consequently give much better 
defined probability functions for the stellar fraction (we 
point out specific examples below). We use the measured 
uncertainties (Table fl]) on the fractions (symmetrized) 
and the total flux to calculate this. Our full expression 



P{ohs,) = B{^\t,-U\\ ]^^ZW (ig) 



.l^k 



C^frac,,fc CrFx,tot,i 



where t/g is the epoch of observation k and P is a nor- 
malization constant such that Pfobs^) = 1. 

Using Equations ([9|, ([l3| and (|T5|, we can express 
Equation (pi) in terms oi Dbservaoles, the microlens- 
ing magnification map histograms, and the intrinsic flux 
probability from the deep field quasar number counts: 



(17) 



These probabilities are plotted for each of the 14 lens- 
ing galaxies in Figure [7| as functions of the stellar mass 
fraction. 

The effect of poorly constrained image fluxes is eas- 
ily seen in the nearly flat probability distribution of 
SDSS 1138+0314, which has only one Chandra observa- 
tion, in which the average uncertainty of the image frac- 
tions is -70%. Compare this to HE 0230-2130, which 
also has only one Chandra observation, but in which the 
average image fraction uncertainty is '^20%. 

5.3. Combined Analysis 

We would like to consider each lensing galaxy as a typ- 
ical member of an ensemble, each with roughly the same 
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Figure 8. Top: distances of quasar images (circles) and 
their mean (stars) from center of lensing galaxy. The red- 
shifts of the lensing galaxies of HE 1113-0641, H 1413+117, and 
WFIJ2026— 4536 have not been measured and were taken to be 
0.7, 0.8, and 0.4, respectively. Their symbols are shown in outline. 
The LM image of RX J0911+0551 is at a radial distance of 17 kpc 
and is not shown. Bottom: distribution of mean radial distance of 
images, Rc, for the 11 lensing galaxies with known redshift. 



configuration such that we are probing the matter con- 
tent at roughly the same radial distance R from the cen- 
ter of the lensing galaxy. To calculate these distances 
in physical units, we use the angular measurements of 
the images and galaxies available on the CASTLES Web 
sitqj along with the redshifts to the lensing galaxies, zi. 

We take the arithmetic mean of the four impact pa- 
rameters where the images form in the lensing galaxy as 
a characteristic radial distance, i?c- Most of the systems 
have Rc within a factor of a few of each other except for 
Q 2237+0305, in which the images form at a mean Rc of 
0.7 kpc, about an order of magnitude less than the mean 
Rc of 6.6 kpc of the other systems (see Figure J8|. We ex- 
clude Q 2237+0305 from the rest of the analysis. We note 
that, had we used the geometric means instead, the num- 
bers would be very similar. The images in Q 2237+0305 
form at a geometric mean radial distance of 0.7 kpc, and 
the geometric mean of the radial distances for the rest of 
the ensemble is 6.1 kpc. 

Unfortunately, zi is not known for HE 1113— 0641, 
H 1413+117, and WFIJ2026-4536 so R cannot be calcu- 
lated for the images of these systems. There have been in- 
dications of lens es of H 1413+117 at redshifts of 0.8, 1.4, 



Kneib et al. 1998 Faure et al. 



and 1.7 (Magai n et al 119881 

|2004( ). ^Morgan et al.| ( |2UU4[ )' estimate a redsh'i'tt o l UTTTor 
the lensinff galaxy or Wl^ l J2026— 4536, and Blackburne 



Tensing g alaxy 
[eFaL] ([2008^ estimate zi = 0.7 for HE 1113-0641. These 
values are all comparable to the redshifts of the other 
lensing galaxies, unlike Q 2237+0305 with zi = 0.04, 
and it is reasonable to assume that their impact param- 
eters are also comparable. We therefore include them 



^ http: //www. cf a. harvard. edu/castles/ 
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Figure 9. Overall probability distribution for the percentage of 
matter in stars including all the X-ray observations for 13 quadru- 
ple lens systems (we do not include Q 2237+0305 — see Section js^ . 
The most likely value for the stellar contribution is 6.8% at a mean 
impact parameter of 6.6 kpc. 
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Figure 10. Normalized probability distributions for the stellar 
fraction in the lensing galaxy of RXJ1131 — 1231 based on split- 
ting the observations into three groups. The first group, shown 
in brown, contains only the observation from 2004 (ObsID 4814). 
The other two groups contain the observations from 2006-2008, 
split into whether the LS image fraction (given in Table [ij was 
higher (blue) or lower (teal) than 0.05. 



in the joint analysis. The radial distances of the im- 
ages in these systems are shown with outlined symbols 
in the top of Figure [S] assuming zi of 0.7, 0.8, and 0.4 
for HE 1113-0641, H 1413+117, and WFIJ2026-4536, 
respectively; they are not included in the histogram in 
the bottom panel of Figure [S] 

For the 10 systems other than Q 2237+0305 with 
known zi , their mean impact parameters Rc are within a 
factor of 2.5 of each other. If we consider all individual R 
in these 10 systems, the spread is nominally a factor of 14. 
Excluding the two extrema (the LS image in B 1422+231 
dit R= 1.2 kpc and the LM image in RX J0911+0551 at 
R = 17 kpc), the spread in R among the 10 systems with 
known zi is only a factor of 3.2. As we discuss in Section 
|6.2[ this is a small enough range in Rc and R that the 
ratio of stellar matter to dark matter is expected to vary 
by only 1.6 over this interval, and we feel comfortable 
combining the individual results to obtain an ensemble 
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Figure 11. As in Figure|9]but also excluding RX J1131-1231. 

result for a mean Rc of (Rc) =6.6 kpc. 

We form the joint probability function of the 
ensemble — excluding Q 2237+0305 — by multiplying to- 
gether the individual probability functions of the 13 lens- 
ing galaxies (shown in Figure ul, and normalizing. The 
results are displayed in Figure^ which shows the joint 
probability distribution of the ensemble for the percent- 
age of matter in stars at a mean impact parameter of 6.6 
kpc. The highest peak of this discrete distribution oc- 
curs at 6.8% stellar matter (93.2% dark matter), and the 
interpolated peak occurs at 6.3% ± 0.3% stellar matter 
(93.7% dark matter). 

6. DISCUSSION 

6.1. RXJ1131-1231 

When exploring the results of the individual observa- 
tions shown in Figure[7[ we noticed that RX J1131 — 1231 
had one of the highest probabilities for a 100% stellar 
fraction, after Q 2237+0305. This is surprising given 
that the first Chandra observation of RXJ1131 — 1231 
displayed a strong signature of significant dark matter 
presence: a hig!:hly suppressed saddle-point image (Black- 1 



burne et al.||2QQ6[ ). We examined the RX J1131-1231' 
probabilities on an observation by observation basis and 
found that, indeed, the first observation (ObsID 4814) 
strongly favored a stellar fraction of 22%. The other 
observations favored higher stellar fractions, either with 
a roughly flat distribution above 22% stars or a strong 
peak at 100% stars. 

We noticed a correlation that those observations with 
a flat distribution above 22% were the ones that had 
lower LS fractions, and the handful of observations (six) 
that peaked at 100% stars were the ones with an LS 
fraction >0.05 in Table [l] We separately analyzed these 
two groups, and the results are shown, along with ObsID 
4814, in Figure [To] 

RXJ1131 — 1231 has the most X-ray observations and 
displays interesting behavior. The HS image evolved 
from being strongly demagnifled by microlensing to be- 
ing strongly magnifled by microlensing, and the LS image 
shows microlensing variations of over a factor of two. It 
may be that our snapshot analysis is not appropriate for 
such complex behavior. Our treatment of each observa- 
tion separately and combination of their weighted results 
discards information on temporal evolution. An analysis 
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that assesses the probabihty of a certain stehar fraction, 
Sj , to produce the entirety of the observat i ons, s imilar to 
the Bayesian III method in Pooley et al. (2009) may be 
more appropriate but is beyond the scope of this work. 

Although we have some minor concerns about the ro- 
bustness of the RXJ1131 — 1231 stellar fraction proba- 
bilities, their effect on our joint analysis shown in Fig- 
ure [9] is minor. If we perform the analysis without 
RX J1131 — 1231, we see that the most probable stellar 
fraction is slightly lower (the interpolated peak is at 
4.6% ± 0.2%), and the probability of 100% stars is near 
zero (Figure 11). 



6.2. Dark Matter Fraction versus Radial Distance 

Our measurements of the stellar mass fraction pertain 
to the impact parameter, that the quad images make 
with respect to the lensing galaxy, and span a range of 
^3-11 kpc. For any given impact parameter, i?, a large 
range of radial distances in three dimensions (i.e., for 
all r > R) is probed within the lensing galaxy. Since 
the percentage of mass in stars is expected to decrease 
with increasing r, we would like to ascertain whether 
our ensemble average likelihood distribution for the dark 
matter fraction (see Figure |9| is well defined, or whether 
we should expect to see a decreasing progression of star 
fraction with increasing mean impa ct para meter, R^. 

Following Koopmans et al.' ("2009) and Schwab et al. 
([2010), we express the three-dimensional light density in 
an elliptical galaxy as 



I{r) = /so 



^0 



(18) 



where /so and ro are constants for a given gala xy, and S 
i s a ni ore nearly universal constant which , Schwab et al. 
(2010) determined to be 

(^ = 2.4 ±0.11 (19) 



based on 5 4 lenses fro m the SLAGS survey (e.g. , Treu| 



et al.||2006l [Koopmans||2006[ [Gavazzi et ar||2007 D. TEe 
value of 0.11 is supposed to represent the rms variation 
in S among different galaxies, rather than an uncertainty 
in the mean value of S. We assume that this power law 
holds over the radial interval r 1-10 kpc. We also 
assume that the stellar mass function and evolutionary 
states of the stars are distance-independent, so that /(r) 
also represe nts the stellar mas s dens ity. 



Similarly, Koopmans et al. (2009) and [Schwab et aT 
( [2010p took the total mass density to be of the form: 



p{r) = Po 



ro 



(20) 



with constants that are analogous to those in Equation 
(4); a is found to be 



a = 1.96 ±0.08 



(21) 



again based on the SLAGS survey (e.g . ,|Treu et al.|2006| 
Koopma ns||2006l [Gavazzi et al.| 2007} koopmans et al. 
2009; Sc hwab et al.||2010p . Similarly, the value of 0.08 is 
supposed to represent an rms variation from galaxy to 
galaxy, rather than an uncertainty in the mean. In this 
expression, p represents both the dark matter and stellar 
contributions to the mass density. 
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Figure 12. As in Figure [9] for systems with known zi, separated 
into two groups based onT^c- Q 2237+0305 and RXJ1131-1231 
have been excluded. 

The observations determine the most probable stellar 
fraction 6', which is the fraction in stars of the total col- 
umn density Ctot along the line of sight at impact pa- 
rameter R. We can integrate expressions (18) and (20) 
to obtain: 

Cstars Jso m-l)/2) T{a/2) /ro^-" 

" PO r((a-i)/2)r(V2) \rJ ^ ' 




where the only radial dependence is in the final factor, 
jl^-s YoT the nominal values for S and a listed above, 
this reduces to 



S{R) 



Cu 



0.77 



ho rro\ 



0.44 



(23) 



Therefore, for a range of mean impact parameters from 
Rc = 3.9 to 9.5 kpc, we expect the stellar fraction S to 
vary by only a factor of ^1.5 due to the dependence on 
the impact parameter. Given that this is roughly the 
resolution of our logarithmic grid of a dozen values of 
S and that our sample is modest in size, we would not 
expect our results to be sensitive to the range in Rc. 

Nevertheless, we divided the 10 lens with known zi 
into two groups: those with Rc < 7 kpc and those 
with Rc ^ 7 kpc. The second group initially contained 
RX J1131 — 1231, but we removed it from t he f ollowing 
analysis given the issues discussed in Section [Oj We ran 
a joint analysis separately on these two groups, and the 
results are shown in Figure [12] As expected, the group 
with smaller (Rc) has a most probable stellar fraction S 
that is larger than the group with larger (Rc)- 

To see how well this result agrees quantitatively with 
Equation (23), we use a fitting function to determine the 
precise location of the peak of the probability distribu- 
tion of each group. The first group, with (Rc) =4.9 kpc, 
peaks at = 5.9% ± 0.5%. The second group, with 
(Rc) = 8.4 kpc, peaks at S2 = 4.5%±0.2%. Based on the 



impact parameter ratio of 1.7, Equation (23) predicts a 



stellar fraction ratio of 1.3. Somewhat reniarkably, given 
the modest size of our samples, S1/S2 = 1.3±0.13. 

We also note the most likely stellar fraction of 100% 
for Q 2237+0305 (see Figure ItI). The lensing galaxy in 
this system is much closer than the others at zi = 0.04 
and consequently has a much smaller Rc of 0.7 kpc. It 
is the only system with a most likely stellar fraction of 
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100%, and this is in quahtative agreement with expec- 
tations. Quantitatively, it is larger than suggested by 



Equation 



parameter. 



but this may be due to either of Equa- 



tions (18) or(20) not being valid at such a small impact 



7. SUMMARY 



We have analyzed 61 publicly available Chandra obser- 
vations of 14 quadruply lensed quasars. We extensively 
tested several methods to reduce and fit the Chandra 
data to obtain the best measurements of the individual 
X-ray fluxes of the q uasar images. As we have shown 



in our previous work ( [Pooley et al.||2007' [2009|), the X- 

oiF 



ray fluxes are a relatively clean measure of microlensing 
effects, unencumbered by source size considerations. 

The results of our data reduction and analysis were 
used in a Bayesian analysis of custom microlensing mag- 
nification maps which marginalized over all observational 
uncertainties as well as multiple observations of a lensed 
quasar. Our analysis yields a most likely local stellar 
fraction of 6.8% (i.e., a most likely dark matter fraction 
of 93.2%) for the ensemble of lensing galaxies, integrated 
along the line of sight at a mean impact parameter of 6.6 
kpc. This is simil ar to the value of 5% found by Medi- 



avilla et al. (2009|), who studied flux ratios in the optical 
and assumed a source size of 2.6 x 10^^ cm. It is als o 
consistent with the recent work of Bate et al. (2011), 
which considered optical and infrared data and found 
dark matter fractions of 50l^S%, 8011°%, and <50% in 
MGJ0414+0534, SDSS J0924+0219, and Q 2237+0305, 
respectively. Those authors performed a marginalization 
over the source size parameters in their analysis. A dis- 
tinct advantage of the work presented here is that our 
X-ray analysis is unencumbered by source-size consider- 
ations. 

We formed two subsets of the lensing galaxies based on 
the mean impact parameters where their images formed 
and found that their most likely stellar fractions varied 
both qualitatively as expected — higher stellar fractions 
closer to the centers of the lensing galaxies — and quan- 
titatively as expected. In addition, we find a most likely 
stellar fraction of 100% for Q 2237+0305, which has a 
mean impact parameter about an order of magnitude 
smaller than all of the other lens systems we studied. 

Our measurement of integrated stellar fraction as a 
function of impact parameter opens up the possibility of 
mapping out the dark matter content of lensing galax- 
ies in a direct and straightforward manner, with minimal 
assumptions, based solely on high quality X-ray observa- 
tions of lensed quasars. 
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